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METHODE POUR MODELISER DES DEPLACEMENTS DE FLU I DES DANS UN MILIEU POREUX. 



(57) - Methode de modeiisation pour optimiser de fa$on 
plus rapide et plus realiste les conditions de deplacement 
dans un milieu poreux rnouillable par un premier fluide (de 
I'eau par exemple), d'un melange de fluides incluant ce flui- 
de mouillant et au moins un autre fluide (de I'huile et even- 
tuellement du gaz). 

- La methode comporte la determination experimental 
de la courbe de variation de la pression capiilaire dans les 
pores en fonction de la saturation relativement aux phases 
liquides. ia modeiisation des pores du milieu poreux par une 
distribution de capillaires a repartition fractale en conside- 
rant, dans le cas d'un melange triphasique eau (fluide 
mouillant) -huile-gaz par exemple une stratification des 
constituants a I'inteheur des pores, avec I'eau au contact 
des parois, le gaz au centre et I'huile formant une couche in- 
tercalate, la determination, a partir de cette courbe de pres- 
sion capiilaire, des valeurs de dimension fractale 
correspondant a une serie de valeurs donnees de la satura- 
tion relativement a la phase liquide, la modeiisation des per- 
meabilites relatives directement sous forme d'expressions 
analytiques dependant des differentes valeurs de dimen- 
sion fractale obtenues et en accord avec la repartition stra- 
tifiee des differents fluides dans les pores, et I'utilisation 
d'un simulateur de milieu poreux pour determiner, a partir 



des permeabilites relatives, les conditions optimales de de- 
placement des fluides dans le milieu poreux. 

- Applications a la production petroliere, a la depollution 
des sols, etc. 
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La presente invention concerne une methode pour modcliser des ecoulements 
di ou triphasiques dans un milieu poreux, basee sur une representation fractale de 
celui-ci. 

La methode selon l'invention est applicable dans de nombreux domaines ou 
Ton veut modeliser les ecoulements de fluides dans des milieux poreux, dans le but 
d'optimiser les conditions de leur deplacement. Parmi les domaines d'application, on 
peut citer par exemple : 

- I'exploitation d'un gisement petrplier et notamment * la production assistee 
d'hydrocarbures par injection de fluides, en utilisant par exemple des injections 
alternees de bouchons de liquide et de gaz (methode dite WAG); 

b) la depollution des sols et notamment des sites industriels, par injection dans 
des couches polluees, de substances telles que des tensio-actifs: 

c) le nettoyage de cakes de filtration par deplacement des substances qui y 
sont retenus; 

d) le sechage du bois; 

- Toptimisation de reactions chimiques par deplacement par exemple de produits 
reactionnels dans une masse de catalyseur, pour augmenter les surfaces de contact, 
etc. 

La methode selon Tinvention est directement applicable par les ingenieurs de 
reservoir pour determiner par exempie la methode de recuperation assistee la plus 
appropriee & appliquer a un gisement souterrain d'hydrocarbures. La methode peut 
egalement servir dans le cadre d'operations de depollution de terrains industriels par 
exemple. 
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Art anterieur 

1) Etudes experimentales: 

La determination experimental des permeabilites relatives d'un milieu poreux 
ou se deplace un fluide polyphasique n'est pas une tache aisee. On simplifie 
5 habituellement les operations de mesure en considerant que Tune des phases est 
immobile dans un etat de saturation irreductible. 

Les valeurs sont acquises par exemple par une methode experimental connue 
dite de "steady state" pour determiner les permeabilites relatives consiste a laisser 
s'ecouler un fluide triphasique avec des rapports de debit imposes entre les phases. 
10 Les permeabilites relatives exprimees en fonction des deux saturations, sont calculees 
en appliquant la loi de Darcy a chaque phase. II est loin d'etre etabli que les mesures 
des permeabilites relatives obtenues par cette methode soient bien representatives des 
deplacements des fluides et, de toute fa?on, elles sont longues a obtenir car a chaque 
changement de regime, il faut attendre qu'un regime d'equilibre setablisse. 

15 Une autre methode connue consiste a realiser des essais en laboratoire pour 

determiner des tableaux de mesures (tel que celui de la Fig. 1 ) reliant les permeabilites 
relatives et les saturations pour chaque couple de fluides du melange triphasique. Par 
ajustement des courbes de production experimentales, on essaie progressivement 
d'ajuster les permeabilites relatives triphasiques. Ces tableaux de donnees sont 

20 introduces alors dans un simulateur numerique de type Score® ou Athos® qui 
calcule les productions de fluides. Cette methode etant basee sur I'acquisition 
prealable de nombreuses mesures experimentales ajustees progressivement par calage, 
est longue a mettre en ceuvre. 

Une autre methode experimental connue de mesure est decrite par exemple 
25 dans les brevets FR-A- 2 708 742 ou 2 724 460 du demandeur (N.B. Aff, 355L 
3693). Elle est mise en ceuvre en pla^ant dans une cellule de mesure un barreau taille 
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dans le milieu poreux a etudier contenant a des saturations variables pour chaque 
experience, deux phases d'un fluide triphasique. On injecte une troisieme phase de ce 
fluide a une extremite du barreau et Ton relcve en fonction du temps la production des 
trois phases a une extremite opposee. 

5 2) Modules 

Le modele empirique connu dit de Stone permet, par des correlations 
empiriques, de predire des donnees relatives a un ecoulement triphasique a partir de 
donnees correspondant a un ecoulement diphasique. II est valide seulement dans les 
cas d'une forte mouillabilite a 1'eau et on le considere generalement comme faiblement 
10 predictif. 

On connait deux types de modeles physiques pour modeliser ies ecoulements 
triphasiques, bases sur Ies courbes de pression capillaire. Celles-ci component 
generalement un pseudo-plateau de faible pente correspondant aux pores les plus gros 
done aux pressions de mercure les plus faibles, et une ou plusieurs parties plus ou 
15 moins lineaires de pente importante. Les courbes de pression capillaire sont en 
relation avec une saturation (par exemple celle du mercure injecte) et un rayon de 
pore, auquel s'arrete le mercure pour une pression d'injection donnee, determinee par 
la loi de Laplace, Pinj.= 2s. r. 

Un premier modele de representation des milieux poreux, est decrit par : 

20 - Burdine, N.T.: "Relative permeability calculations from pore size distribution data" 
Trans AIME (1953) Vol. 198; ou par 

- Corey, A.T.: "The interrelation between oil and gas relative permeabilities". Prod 
Monthly (1954) Vol. 19,38. 

Suivant ce modele, le milieu poreux est represente par un faisceau de 
25 capillaires cylindriques avec une distribution de rayons qui sont donnes par la courbe 
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de pression capillaire obtenue par injection de mercure. Lcs permeabilites sont 
obtenues par application de la loi de Poiseuille a lecoulement do fluides dans ces 
capi II aires. 

Ce modele est base sur la representation du milieu poreux comme un 
5 assemblage de capillaires dc rayons differents. La relation entre le volume et le rayon 
des pores est donnee par la valeur de la pente du pseudo-plateau. Les trois fluides 
sont censes se partager les capillaires, le fluide mouillam (l'eau) occupant les plus 
petits, le fluide le moins mouillant (le gaz), les plus gros. le troisieme fluide (1'huile) 
une zone de pores de dimensions intermediaires. II n'est pas possible de decrire les 
10 interactions entre les fluides puisque dans un tel modele, ils secoulent dans des 
chenaux separes. Ce modele enfin, ne peut etre utile que si lc pseudo-plateau couvre 
une large gamme de saturations. Suivant ce modele, les trois phases d'un ecoulement 
triphasique se deplacent dans des capillaires differents et il n'y a entre elles aucune 
interaction. 

15 Un autre modele physique connu de representation d'un milieu poreux est 

decrit par : 

- de Gennes, P.G.: "Partial Filling of a Fractal Structure by a Wetting Fluid" Physics 
of Disordered Materials 227-241, New York Plenum Pub. Corp. (1985); ou par 

- Lenormand, R. "Gravity Assisted Inert Gas Injection: Micromodel Experiments 
20 and Model based on Fractal Roughness", The European Oil and Gas Conference 

Altavilla Milica, Palermo, Sicily, (1990). 

Suivant ce modele, on considere que la surface interne des pores est isotrope 
et a caractere fractal et peut etre modelisce sous la forme d'un "bouquet" de rainures 
capillaires paralleles qui font que les pores presentent une section transversale fractale. 
25 La section transversale de chaque pore est construite suivant un processus iteratif 
(Fig.l). On divise le demi-perimetre d'un cercle de rayon R n en T| parties et en 
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remplafant chacune de ces r| parties par un demi-cercle ou rainure. A chaque etape k 
du processus, N k nouvelles rainures semi-circulaircs dc rayon R k et de section totale 
Aw sont creees. 

La dimension fractale DL de la section transversale a ia fin de I'etape k est 
reliee au nombre d'objets Nk genere avec l'echelle donnee I K , par la relation : 

La dimension fractale peut etre deduite dune courbe de pression capillaire de 
mercure selon la procedure suivante. Du mercure est injecte dans un milieu poreux 
avec une pression d'injection qui croTt par paliers. La loi de Laplace permet de deduire 
le volume des pores connaissant le volume de mercure injecte pour une pression 
d'injection donnee et Ton peut construire la courbe de pression capillaire en drainage 
reliant la pression d'injection a la quantite de mercure injectee et la courbe reliant la 
proportion du volume total occupe par les pores et la taille des pores. Dans le cas ou 
Ton draine un liquide mouiliant hors du milieu poreux tel que de I'eau par une injection 
de gaz, la correlation entre la pression capillaire gaz-eau et la saturation de la phase 
mouillante est donnee par : 

i 

Dans cette representation log-log, la pente l/(D L -2) de chaque ligne (de -1.5 a 
-3.3) donne la valeur de la dimension fractale lineaire D L entre 1 ,3 et I J. 

Les resultats experimentaux montrent aisement que les valeurs des 
permeabilites relatives gaz-eau exprimees en fonction des trois saturations, obtenues a 
partir des expressions qu'en donnent les modeles connus et les modes de repartition 
des phases dans la structure des pores, sont loin des valeurs mesurees et done que les 
modeles en cause se revelent par trop simplistes pour representer les interactions 
complexes qui se produisent entre les phases fluides. 
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Definition de la methode selon I'invention: 

La methode de moderation selon I'invention permet d'optimiser de fa<;on plus 
rapide et plus realiste les conditions de deplacement dans un milieu poreux mouillable 
par un premier fluide, de melanges di ou triphasiques incluant ce premier fluide 
mouillant et au moins un deuxieme fluide. Elle offre de ce fait aux operateurs, un outil 
plus fiable pour evaluer notamment les meilleurs modes de deplacement des fluides au 
sein du milieu poreux. 

Elle est caracterisee en ce qu'elle comporte en combinaison : 

- la determination experimentale de la courbe de variation de la pression capillaire 
(P c ) dans les pores d'un echantillon de ce milieu poreux en presence dun fluide 
mouillant et d'un fluide non mouillant, (par injection de mercure dans un 
echantillon place sous vide par exemple); 

- la moderation des pores du milieu poreux par une distribution de capillaires a 
section fractale en considerant une repartition stratifiee des fluides a l'interieur des 
pores, le fluide mouillant se repartissant au contact des pare is et autour du 
deuxieme fluide (ou des deux autres en triphasique); 

- la determination a partir de cette courbe de pression capillaire, des valeurs de la 
dimension fractale correspondant a une serie de valeurs donnees de la saturation 
relativement aux liquides; 

- la modelisation des permeabilites relatives directement sous forme d'expressions 
analytiques dependant des differentes valeurs de dimension fractale obtenues; et 

- Introduction des permeabilites relatives dans un simulateur du milieu poreux et la 
determination au moyen de ce simulateur, des conditions optimales de deplacement 
des fluides du melange dans le milieu poreux. 




m 

7 



La methode s'applique par exemple a, la determination des deplacements de 
melanges de fluides comportant un premier tluide mouillant, un deuxieme fluide non 
mouillant et un gaz, en considerant une repartition stratifiee des tluidcs a 1'interieur 
des pores, le fluide mouillant se repartissant au contact des parois, le gaz occupant le 
5 centre des pores et le deuxieme fluide etant reparti sous la forme dune couche 
annulaire au contact a la fois du gaz et du premier tluide. 

La methode peut s'appliquer notamment a la determination par un simulateur 
de gisement, des caracteristiques optimales de substances additionnees a des 
bouchons de fluide mouillant injectes dans une formation en alternance avec des 
10 bouchons de gaz, dans le but de deplacer des hydrocarbures en place, ou bien encore 
celles dun fluide injecte dans le sol pour deplacer des substances polluantes. 

La modelisation des phenomenes obtenue par la methode presente de grands 
avantages. Elle permet une meilleure adequation avec les resuitats obtenus en 
laboratoire car les phenomenes physiques sont mieux pris en compte. Les resuitats du 
15 modele sont meilleurs de ce fait lors d'un changement dechelle par exemple, pour la 
modelisation d'une application sur un champ d'operations. 

Le temps de calcul est reduit compare a celui qui est necessaire quand on 
utilise des tables, comme avec les methodes anterieures. La modelisation de type 
fractal realise est mieux a mcme de traiter les effets d'hysteresis rencontre dans 
20 ('utilisation des procedes d'injection de type WAG. 

Les resuitats de la methode sont en outre parfaitement integrables dans de 
nombreux simulateurs de reservoir : simulateurs 3D, heterogenes, compositionnels, 
etc. 

Sexploitation des resuitats par les logiciels duplication est facilitee. II n'est 
25 plus necessaire de faire des interpolations risquees comme e'est la regie quand on 
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opere a partir des valeurs discretes des tableaux de resultats, pour tracer des isoperrns 
par exemple. 

D'autres caracteristiques et avantages de la methode selon I'invention, 
apparaitront a la lecture de la description ci-apres dun exemple non limitatif de 
5 realisation, en se referant aux dessins annexes ou : 

- la Fig.1 qui illustre sous forme dun tableau, les liens existam pour un melange 
triphasique entre les valeurs obtenues experimentalement. de permeabilites 
relatives d'un fluide et les saturations pour deux des trois fluides; 

- la Fig.2 montre une representation fractale d'un pore; 

la Fig.3 montre de fa^on schematique la repartition des phases d'un fluide 
triphasique dans un pore fractal avec le fluide mouillant W au contact de la paroi, 
la phase gazeuse G repartie dans la plus grande partie du volume du pore (rayon 
R G ), le fluide non mouillant O etant une couche entre le fluide mouillant et le gaz; 

la Fig.4 montre en fonction de la saturation en eau, la courbe de pression 
capillaire d f un echantillon de gres des Vosges par exemple, dont la pente locale S 
sert a determiner la dimension fractale des pores; et 

- la Fig.5 montre en fonction du temps, les courbes de production obtenues 
experimentalement pour du gaz (G) pour de l'eau (W) et pour de I'huile (O), 
comparees aux courbes equivalentes obtenues par simulation en utilisant la 

20 methode selon 1'invention. 

Description detail lee de la methode 

La methode selon Tinvention permet de determiner les permeabilites relatives 
triphasiques de milieux poreux, en s'appuyant sur un modele de type fractal du milieu 
poreux, en se basant sur une approche decrite par : 



10 
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- Kalaydjian, F.J-M et al. : "Three Phase Flow in Water-wet Porous Media: 
Determination of Gas-oil Relative Permeabilities under Various Spreading 
conditions"; 68 lh Ann. Tech. Conf. And Exh. of the SPE. Houston, Texas, 1993. 

Description de l'invention 

La methode selon l'invention comporte, comme on I'a vu une modelisation de 
I'ecoulement des phases avec repartition des phases a I'interieur de la structure fractale 
de chaque pore. Dans le cas d'un ecoulement d'eau et d'huile au sein d'une structure 
po reuse mouillable a I'eau, I'J^ pore. Dans le cas d'un 

ecoulement triphasique, il y a repartition stratifiee, I'eau qui est le fluide mouillant, 
s'ecoule le long des parois des pores, le gaz circule dans ie volume du pore et I'huile 
s'ecoule entre le gaz et I'eau. Les saturations sont calculees comme etant la surface 
relative dans une section trans versale occupee par chacun des fluides. 

A l'equilibre, toutes les rainures ayant un rayon inferieur ou egal a Rk qui est 
donne par la loi de Laplace Pc=2 y/R k , sont occupes par le fluide mouillant, les tubes 
les plus gros, par le fluide non mouillant. Ainsi on exprime la saturation en fluide 
mouillant comme la fraction de I'aire des tubes occupes. 

En calculant la fraction de I'aire des capillaires occupee par I'eau pour tous les 
rayons entre R k et Rpo , on aboutit a l'expression : 




et comme Pc=2 y/R k , la correlation entre la pression capillaire et la saturation de la 
phase mouillante est donnee par 

ou Sw est la saturation de la phase mouillante. 



2772483 



10 

La representation graphique de cette correlation dans un diagramme log-log 
est une ligne droite partant du point (Sw, Pc) correspondant au capillaire le plus grand 
de la structure fractale avec un rayon R 0 . On peut supposer que (Fig.4), 

le rayon Ro, le premier envahi quand le mercure est injecte, correspond a une 
valeur de la saturation de l'ordre de l/r=10 3 . Chaque segment de la courbe de 
pression capillaire est une partie d'une ligne partant de R 0 (supposee la meme pour 
tous les differents segments), correspondant a la correlation P c , S w deja cit£e. 
Chaque ligne a une pente donnee, on peut lui associer une dimension lineaire 
fractale. Les valeurs de la pente varie de -1,5 a -3.3 comme le montre la Fig.4, ce 
qui conduit a des valeurs de la dimension lineaire fractale D L entre L3 et 1.7. 

- chaque domaine est atteint par le mercure pour des saturations correspondant a 
l'emplacement ou on trouve Ro sur chaque ligne. 

La saturation des deux liquides quand le gaz est present dans le pore est 
calculee comme expliquee plus haut pour une phase : 

at 

en supposant que les deux liquides occupent les tubes dont les rayons sont inferieurs 
ou egaux a R K et le gaz, le centre de chaque pore. La saturation en huile est Taire 
relative de la section transversale des rainures occupee par I'huile dont le rayon est 
inferieur ou egal a R K . 

20 Pour proceder au calcul des permeabilites relatives, il est necessaire de 

determiner pour chaque phase la fraction qui circule et par consequent, d'estimer de 
facpon systematique, les saturations correspondant aux fractions stagnantes. Ceci doit 
etre fait pour les deux cas etudies, par exemple le drainage de Teau et de Thuile pur le 
gaz, et rimbibition en eau. 
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Calcul des permeabilites relatives: 

Permeabilites relatives des liquides 

L'application de la loi de Poiseuille a chaque capillairc du faisceau occupe par 
la phase qui I'occupe, permet de calculer les permeabilites relatives a I'eau et Thuile 
(Fig.5). 

En ne considerant que la fraction en circulation qui contribue a la conductivity 
hydraulique, les permeabilites relatives pour I'eau et I'huile s'expriment de la maniere 
suivante : 

k„=sr*-sf^ (i) 

= k nt (2r».)[ - (s w + s or y^ ] (2) 

Dans ces expressions, il est utile de mentionner que 

- la saturation en eau irreductible Swi est supposee stable. 

- Sorw represents la saturation en huile residuelle laissee en place apres un balayage 
du milieu a I'eau et Sor est une partie de cette saturation residuelle qui correspond 
a une saturation en eau determinee. 

- La gamme de tailles de capillaires occupee par Thuile mobile est calculee comme la 
difference entre les tailles des capillaires occupes par les deux liquides avec la 
saturation totale en liquide S L =S 0 4-S W , et celles des capillaires satures en I'eau et en 
huile stagnante. 

- Kro(2ph.) represente la valeur de la permeabilite relative a i'huile determinee par un 
test d'imbibition a I'eau et I'huile. Quand I'eau et Thuile sont seuls presentes et 
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comme l'echantillon teste est mouillable a I'eau, I'huile secoulera dans la section du 
pore exactement comme le fait le gaz dans un ecoulement triphasique. 

Permeabilite relative au gaz 

Comme ie gaz est une phase non mouillante, il occupe I'espace central du pore 
5 et il s'etend vers sa peripheric au fur et a mesure que la saturation en gaz augmente, 
sans toutefois arriver au contact de la paroi solide (Fig. 6). On considere que le gaz 
circule dans un seul pore dont le rayon R g est donne par la relation : 

R g = Ro+R,+R 2 + ... + R k , 

la permeabilite au gaz etant donnee alors par : 

10 * nr =^«.(l-5?) 4 (4) 

avec a = 2^7 * Pl representant la dimension fractale lineaire du milieu poreux et S L 
* etant la saturation totale en liquide egale a 1-S g . 

Ayant determine les permeabilites relatives a Teau, a I'huile et au gaz, 
respectivement k^, k r0 et krg, on utilise un simulateur de reservoir adapte aux 
15 conditions de laboratoire tel que le simulateur ZCORE®, par exemple, suivant une 
procedure connue que Ton trouve decrite par exemple par : 

- Moulu J.C. et al: "Performance and Numerical Interpretation of Gas Drainage 
Core Tests under Secondary and Tertiary Conditions", Paper SCA 9508, SCA 
Symposium, San Fransisco, Sept. 1995. 

20 On peut done facilement optimiser les conditions a remplir pour deplacer des 

fluides petroliers en place dans un gisement, que ce soit par injection de gaz soit par 
injection alternee de bouchons d'eau et de gaz (methode ditc WAG), en tenant compte 
des conditions de pression et de temperature regnant a la profondeur de production. 



BNSDOCID: <FR 2772483A1_I_> 




2772483 



Validation : 

La methode selon Tinvention a ete validee en realisant des injections de gaz 
dans des milieux poreux contenant de I'eau et de I'huile. dans des conditions tres 
variees. On voit sur la Fig. 5 par exemple que Ton obtient un tres bon accord entre les 
courbes de production des trois phases (eau, huile, gaz) obtenues experimentalement 
et celles predites par le simulateur de gisement alimente par les donnees obtenues en 
accord avec la methode. 
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REVENDICATIONS 

1) Methode de modeiisation pour optimiser les conditions de deplacement 
dans un milieu poreux mouillable par un premier fluide, dun melange de fluides di ou 
triphasiques incluant ce premier fluide mouillant et au moins un autre fluide, 
caracterisee en ce qu'elle comporte en combinaison : 

- la determination experimental de la courbe de variation de la pression capillaire 
(P c ) dans les pores dun echantillon de ce milieu poreux eh presence d'un fluide 
mouillant et d'un fluide non mouillant, 

- la modeiisation des pores du milieu poreux par une distribution de capillaires a 
section fractale en considerant une repartition stratifiee des fluides a I'interieur des 
pores, le fluide mouillant se repartissant au contact des parois et autour du dit au 
moins un autre fluide; 

- la determination a partir de la dite courbe de pression capillaire, des valeurs de 
dimension fractale correspondant a une serie de valeurs donnees de la saturation 
relativement aux phases liquides; 

- la modeiisation des permeabilites relatives directement sous forme d'expressions 
analytiques dependant des differentes valeurs de dimension fractale obtenues; et 

- Introduction des permeabilites relatives dans un simulateur du milieu poreux et la 
determination au moyen de ce simulateur, des conditions optimales de deplacement 
du melange de fluides dans le milieu poreux. 

2) Methode selon la revendication 1, appliquee a des melanges de fluides 
comportant un premier fluide mouillant, un deuxieme fluide non mouillant et un gaz, 
caracterisee en ce Ton modelise les pores du milieu poreux par une distribution de 
capillaires a repartition fractale en considerant une repartition stratifiee des fluides a 
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Tinterieur des pores, ie fluide mouillant se repartissant au contact des parois, le gaz 
occupant le centre des pores et le deuxieme fluide etant reparti sous la forme d'un film 
annuiaire au contact a la fois du gaz et du premier fluide. 

3) Methode selon la revendication 1 ou 2, caracterisee en ce que Ton met en 
5 oeuvre le simulateur de gisement pour determiner les caracteristiques opttmales de 

substances additionnees a des bouchons de fluide mouillant injectes dans une 
formation en alternance avec des bouchons de gaz dans le but de deplacer des 
hydrocarbures en place. 

4) Methode selon Tune des revendications precedentes. caracterisee en ce 
10 qu'elle comporte la mise en ceuvre d'un simulateur de gisement pour determiner les 

caracteristiques optimales d'un fluide injecte dans le sol pour drainer des substances 
polluantes. 
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